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Abstract 

The growth rate of matter perturbation and the expansion rate of the Universe can be used 
to distinguish modified gravity and dark energy models in explaining cosmic acceleration. We 
explore here the inclusion of spatial curvature into the growth factor. We expand previous re- 
sults using the approximation and then suggest a new form, fa = ^Tn + (7 — 4/7)r2fc, as 
an approximation for the growth factor when the curvature il.^ is not negligible, and where the 
growth index 7 is usually model dependent. The expression recovers the standard results for the 
curved and flat ACDM and Dvali-Gabadadze-Porrati models. Using the best fit values of J7m0 
and r^fco to the expansion/distance measurements from Type la supernovae, baryon acoustic os- 
cillation, WMAP5, and H{z) data, we fit the growth index parameter to current growth factor 
data and obtain 7a {^k / 0) = OmtoH and 'Jdgp i^k / 0) = 0.53l:[J:}^. For the ACDM model, 
the 1-a observational bounds are found consistent with theoretical value, unlike the case for the 
Dvali-Gabadadze-Porrati model. We also find that the current data we used is not enough to put 
significant constraints when the 3 parameters in are fit simultaneously. Importantly, we find 
that, in the presence of curvature, the analytical expression proposed for fa provides a better fit 
to the growth factor than other forms and should be useful for future high precision missions and 
studies. 
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I. INTRODUCTION 



The discovery of late time cosmic acceleration [l| challenges our understanding of the 
standard models of gravity and particle physics. Within the framework of Friedmann- 
Robertson- Walker cosmology, a dark energy component with negative pressure is invoked 
to explain the observed accelerated expansion of the Universe. One simple candidate of 
dark energy which is consistent with current observations is the cosmo logical constant. By 
choosing a suitable equation of state w = p/ p for dark energy, we can recover the observed 
expansion rate H{z) and the luminosity distance redshift relation dilz). Many parametric 
and nonparametric model-independent methods were proposed to study the property of dark 
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energy, see for example 
therein. 

The apparent acceleration of the Universe may be explained by the modification of grav- 
itation, such as for example the Dvali-Gabadadze-Porrati (DGP) brane- world model jl^, 
in which gravity appears four dimensional at short distances but is modified at large dis- 
tances. Recently, it was shown that the late cosmic acceleration of the Universe can also be 
realized in the Horava-Witten heterotic M theory {20I and string theory 21] on 5*^/^2. It is 
remarkable that the acceleration is transient in all these models. To distinguish the effect of 
modified gravity from dark energy, we can use the growth rate of large scale in the Universe 
in addition to the distance data. While different models give the same late time accelerated 
expansion, the growth of matter perturbation they produce differs [22|]. Recently, the use of 



the growth rate of matter perturbation in addition to the expansion history of the Universe 



to differentiate dark energy mode 
example, an incomplete list 
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s and modified gravity attracted much attention, see, for 
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and references therein. To the linear order of perturbation, at large scales, the matter density 
perturbation 6 = Spm/Pm satisfies the simple equation: 



6 + 2H6 - AnGes pm5 = 0, 



(1) 



where pm is the matter energy density and Gefr denotes the effect of modified gravity. For 
example, G^s/G = (4 + 2uj)/{3 + 2uj) for the Brans-Dicke theory [41] and G^s/G = 1 + [3 - 
6rcH{l + H /3H^)]^^ for the DGP model j42, 43], the dimensionless matter energy density 



STrGprn/i^H"^)- In terms of the growth factor / = din 6 /din a, the matter density 
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perturbation Eq. ([T]) becomes 

f' + f + 



H 



+ 2 / 



3 Geff ^ 

2 G ™' 



(2) 



where /' = df/dlna. It is very interes ting that the solution of the equation can be ap 



proximated as / = fi^ 
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491 ] and the growth index 7 can be obtained 



for some general models. The approximation wa s first proposed by Peebles for the mat- 
ter dominated universe as f{z 



0) 



44l |: then a more accurate approximation, 



/(^ = 0) 



^Q^'^ , for the same model was derived in 



45I . |46| | . For a dynamical dark energy 



mode 
in 



47 



with slowly varying w and zero curvature, the approximation f{z) = Q{zy was given 



48l | . For more general dynamical dark energy models in fiat space, it was found that 
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0.55 + 0.05[1 + w{z = 1)] with w > -1 and 7 = 0.55 + 0.02[1 + w{z = 1)] with w < -1 



5l|. For the fiat DGP model, 7 = 11/16 |50|. 



Fitting f{z) = Q{zy to several sets of the most recently observational data, lately one of 



us found that ' 
Porrati model 



0.64+^:{^ for the ACDM model and 7 = 0.55l^:t| for the Dvali-Gabadadze- 
491]. In this paper, we shall generalize such studies to the case where the space 



is not fiat. In particular, in Sec. II we consider the dark energy model with constant w, 
while in Sec. Ill we discuss the DGP model. In Sec. IV, we apply the union compilation 
of type la supernovae (SNe) data js^], the baryon acoustic oscillation (BAO) measurement 



from the Sloan Digital Sky Survey 



531 ] ■ the shift parameter, the acoustic scale I a and the 



redshift 2;* of the last scattering surface measured from the Wilkinson Microwave Anisotropy 



Probe 5 yr data (WMAP5) |54| . and the Hubble parameter data H{z) [55|, l56|. We also use 




64| to find the constraint on 



the growth factor data f(z) [38 



the growth index 7. We also fit our simple analytical expression fa = fi^ + (7 — 4/7)fijt to 
all the data sets by allowing the three parameters flmo, ^ko and 7 all to vary. We conclude 
the paper in Sec. V. 

It should be noted that when the Universe is not fiat, the approximation f{z) = + 
^a/SO was proposed in j65] and f{z = 0) = r^^^ + ^]Ao(l + ^^,no/2)/70 in [66] for the ACDM 
model. From these works one can see that the approximation for the growth factor may not 
be simply when the curvature is not zero, and that the growth index 7 can still be used 
as the signature of modified gravity and dark energy models. 
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II. DARK ENERGY MODELS WITH CURVATURE 



For the curved dark energy model with a constant equation of state w, we have 



H 



-a 



2 2 

The energy conservation equation tells us that 



[l+w(l-n^-fifc)]. 



(3) 



(4) 



Substituting Eqs. ([3D and (g]) into Eq. ([2D, we get 



Plugging / = into Eq. ([5D, we get 



11 3 



f = ^^m. (5) 



[3w{l-nm-nk)-^k]^n.lnnm^+ (^7 " ^ [3^(1 - -l^fe) " l^fc] +fi7n " ^ + ^ = 0- 

(6) 

At high redshift, both 1 — Qrn and Qk are small quantities and they evolve differently. If we 
expand Eq. ([6]) around = 1, then we also need to deal with flk ~ If ^k = 0, to 

the first order of (1 — flm), we get [48|, l67l | 

3 



_ (1-^) 3 {l-w){l-3w/2) 
^ 5-Qw 125(1 -6w;/5)2(l - 12w/5)^ 

If we take flk as a constant, then to the lowest order, we get 



(7-l/2)(3ti; + l)(]fc = 0, 



(7) 



(8) 



so 7 = 1/2 which is not consistent with other results. Therefore, we cannot take as a 
constant. If we think that Qk and 1 — Qm are independently small quantities, then to the 
first order of 1 — Qm and Qk, we get 

3{w-l){l-nm) - (3W+ l)fifc 



7b 



(9) 



{6w - 5)(1 - n^) - 2{3w + l)fifc ■ 

If = 0, Eq. (Q recovers the result given by Eq. ([7]) for the flat space. If there is no dark 
energy, = 1 — flm, then Eq. tells us that 7 = 4/7 which is consistent with the results 
obtained in For the ACDM model, we get 7 = [6(1 -2^]fe]/[ll(l -fi^) -40fc]. 

Equation is also consistent with Eq. (15) in 68| giving a very similar expression when 



expanded in terms of Qde = 1 —^m — ^k and Qk- Although the result in Eq. ([9]) is consistent 
with other results, the expression for 7 is perhaps too complicated to be a constant; and 
we derive below a more practical expression for 7 to distinguish between different models. 
This means that the approximation needs to be modified so that we can get a simple 
number to distinguish different models. Let us recall the approximation used for the curved 
model with a cosmological constant. Martel used the approximation / = Q'^ + i^A/SO |65 |. 
In [66|, the approximation / = Q'^ + + fim/2)/70 was used. Following these ideas, we 
propose the approximation fa = fij^^ + Pflk- Now Qm and Qk are independent variables, so 

^7 



(10) 



Substituting the expression of fa, Eqs. (jlj) and ( ITOl) into Eq. ([2]), expanding all quantities 
around Qm = 1 and Qk = 0, and then keeping terms linear in 1 — Qm — and Qk, we finally 
find that 

3(1-^) . 4 
Therefore, the approximation becomes 

fa = ^l + {i-^mnk, 7 = ^^^- (12) 

For this approximation, we have only one constant 7 and this growth index can be used as the 
fingerprint of the model. When = 0, the approximation (IT^ recovers the familiar result. 
For the ACDM model, the approximation is fa = ^ni^^ — 2^7^/77. This result is consistent 



with Eq. (15) in [68| for 7 when the latter is used into their function /(a) = ilrnio-)"' and 
expanded around Qk = 0. 

In order to see how well the approximation fa fits the growth factor /, we need to solve 
Eq. ([2]) numerically with the expression of Qm and Qk- The dimensionless matter density is 

fi^o + fi,o(l + ^)-i + (l-fi^o_nfco) (1 + ^)3-' ^ ^ 

and the dimensionless curvature density is 

" n^o + fi,o(l + z)-l + (l-l]„0-l],o) (1 + ^)3-- ^ ^ 

By using the above Eqs. (IT5]) and ffT^ . we solve Eq. ([2]) numerically to get the growth factor 
/ for different values of ^2^0, ^fco and w. Then we compare the approximation /„ with / 
by plotting the relative error [fa — f)/f for the ACDM model in Fig. [H From there we see 



that the error of the approximation is a few percent, so fa approximates / very welL We 
also compare the approximation Q'J^ with '-fb given by Eq. (Q with the approximation fa by 
plotting the absolute value of the ratio \ ~ f)/{fa ~ f)\ in FigEl from which we see that 
fa approximates / better. For the dark energy model with constant w, the accuracy of the 
approximation fa is shown in Fig. [3] and the comparison between the approximation of /„ 
and fb is shown in Fig. IHfor some typical values of Qmo, ^ko and w. The approximation fa 
is not only a better approximation, but also has a simpler expression for the growth index 
7. Therefore, we should use Eq. (fT2|) to approximate the growth factor. 
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FIG. 1: The relative difference between the growth factor / and the approximation fa with 7 = 6/11 
in Eq. ^ for the ACDM model. 



III. DGP MODELS WITH CURVATURE 



For the DGP model, we have 

Geff ^ ml - 4(1 - + 2v/r^T4(3 - 41],. + + ^1) 

G 31]^-3(l-l]fc)2 + 2v^r^T4(3-4fifc + 2fi^fifc + fi2)- ^ 
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FIG. 2: The ratio \ {0,m — f)/ifa ~ f)\ of the two different approximations of the growth factor for 
the ACDM model. Since the case with ^ImO = 0.3 and ^l^o = —0.1 was cut off in the plot, we show 
it explicitly in the inset. 



The Friedmann equation gives 



(16) 



The energy conservation equation tells us that 



= -fi™ 3 - 2fi 



(17) 



1 + — , 

As discussed in the previous section, we consider the approximation = + Substi- 
tuting the expression of /a, Eqs. ( TT7|1 and (ITUl) into Eq. ([2]), expanding all quantities around 
^rn = 1 and Qk = 0, and then keeping terms linear in 1 — Qm — ^k and Qk, we obtain 



2 + ^(/3-7) 



-47 + - 



, 1 



0. 



(18) 



Therefore, the approximation of the growth factor for the DGP model is 



A - + - nk 



(19) 
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FIG. 3: The relative difference between tlie growth factor / and the approximation fa in Eq. ([12 
for the dark energy model with constant w. 
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FIG. 4: The ratio \ {fa — f)/{^m — f)\ of the two different approximations of the growth factor for 
the dark energy model with constant w. 
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On the other hand, if we approximate the growth factor / by fij^, then set /3 = in Eq. 
(fT8|). finally we find 

^ ll(l-^]^)-3^]fc 
l6(l-^]„)-2^],• ^ ^ 

In the DGP model, the dimensionless matter energy density is given by 



where fi^o = (1 ~ ^mo ~ ^A;o)/2-\/1 ^ ^fco- The dimensionless curvature energy density is 

= + (22) 



Combining Eqs. ([2]), f[T5|) . f[T6|) . fl2T|) and fl22l) . we get the evolution of the growth factor /, 
and the result is then compared with the approximation /„. In Fig. [5l we show the accuracy 
of the approximation by plotting {fa — f)/ f- 




FIG. 5: The relative difference between the growth factor / and fa with 7 = 11/16 for the DGP 
model. 

From Fig. [5l we see that the error is under 10%. We also checked the accuracy of the 
approximation fi^ and find that the approximation fa is usually better for the model we 
are interested in. 
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IV. OBSERVATIONAL CONSTRAINTS 



Now we use the observational data to fit the dark energy model with constant w and 
the DGP model. The parameters in the models are determined by minimizing = xL + 
Xbao + Xcmb Type la SNe data, we use the reduced union compilation of 307 

Type la SNe 52|. The union compilation has 57 nearby Type la SNe and 250 high-^ Type 



69| and the ESSENCE Survey [TO, IZll, 



la SNe. It includes the Supernova Legacy Survey 
the older observed SNe data, and the extended data set of distant SNe observed with the 
Hubble space telescope. To fit the Type la SNe data, we define 

307 



i=l 



a} 



(23) 



where the extinction-corrected distance modulus n^z) = 5\ogiQ[dL{z) /Mpc] + 25, fiobs is the 
observed distance modulus, is the total uncertainty in the SNe data, and the luminosity 
distance is 



1 + z 



:smn 



dz' 



where 



smn 



X) 



sm{^/\Qk\x) 



E z' 



if Qk < 0, 
if = 0, 
if ilk > 0, 



(24) 



(25) 



^smh.{^/\Qk\x) 

and the dimensionless Hubble parameter E{z) = H{z)/Hq = [Qo{l + z)^ + (1 — fio)(l + 
2,-)3(i+«))ji/2 ^Yie dark energy model with constant w and E(z) = [i^o(l + z)^ + (1 — 
fio)V4]i/2 + (1 - Qo)/2 for the DGP model. 



To use the BAO measurement from the Sloan Digital Sky Survey data, we define 531] 



X,bao 



r,{zd)/Dv{z = 0.2) -0.198 



0.0058 

where the effective distance is 



r,{zd)/Dv{z = 0.35) - 0.1094Y 
0.0033 ) 



(26) 



Dv(z) 



dl{z) 



[l + zYHiz) 



1/3 



(27) 



The redshift Zd is fitted with the formulas 721] 

1291(fi„/i2)0-25i 



Zd 



0.659(l]^/l2)0.828 



[l + 6i(fife/i 



2^b2^ 



(28) 
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6i = 0.313((]„/i2)-o-4i9[l + 0m7{nmhy-^'% 62 = 0.238{nrahY^^\ (29) 
and the comoving sound horizon is 



where the sound speed c,(z) = 1/V3[1 + Rb/{1 + z)], Rb = 315000nbh'^{Tcmb/2.7K)-^. 

To implement the WMAP5 data, we need to add three fitting parameters R, la and z^,, 
so xlmb — ^XiCov~'^{xi, Xj)Axj, where Xj = {R, la, z^) denotes the three parameters for 
the WMAP5 data, Axi = Xi — xf"^ and Cov{xi,Xj) is the covariance matrix for the three 
parameters js^. The acoustic scale I a is 



_ TTdiXz^) 



where the redshift z^, is given by 



73| 



z^ = 10A8[1 + 0.0012A{Qbh^)~'^-'^^^][l + gii^mhy^] = 1090.04 ± 0.93, (32) 
_ 0.0783(1] j,/i^)-°-^3g _ 0.560 

~ l + 39.5(a/l2)0.763' 92 - l + 21.1(l]fe/l2)l-81- ^ ^ 

The shift parameter 

R = ^^sinn f ^1 = 1-710 ± 0.019. (34) 

Simon, Verde, and Jimenez obtained the Hubble parameter H{z) at nine different red- 
shifts from the differential ages of passively evolving galaxies 55|]. Recently, the authors in 
3 obtained H{z = 0.24) = 83.2 ± 2.1 and H{z = 0.43) = 90.3 ± 2.5 by taking the BAO 
scale as a standard ruler in the radial direction. To use these 11 H{z) data, we define 

2 _ [Hobs{zi) - H{zi)]'^ 

Xh - -2 ' y'^'^) 

i=l ^hi 

where au is the la uncertainty in the H{z) data. We also add the prior ifo = 72 ± 8 



km/s/Mpc given by Freedman et al. 7^. The likelihood for the parameters in the model 
and the nuisance parameters VLbfi^ and Hq is computed using a Monte Carlo Markov Chain 
(MCMC). The MCMC method randomly chooses values for the above parameters, evaluates 
X^, and determines whether to accept or reject the set of parameters using the Metropolis- 
Hastings algorithm. The set of parameters that is accepted to the chain forms a new starting 
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FIG. 6: The la and 2a contours of 0, and w by fitting the dark energy model with constant w to 
the combined data. The point with -|- denotes the flat ACDM point 0,^ = and w = —1. 
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FIG. 7: The la and 2a contours of and 17^ for the ACDM model. 
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FIG. 8: The la and 2a contours of Qrn and for the DGP model. 



point for the next process, and the process is repeated for a sufficient number of steps until 
the required convergence is reached. Our MCMC code is based on the pubhcly available 
package COSMOMC [TS]. 

By fitting the dark energy model with constant w to the combined data (except the 
growth data that we use next), we get = 332.6, fi^o = 0.276+o:o26, fifco = 0.002 ± 0.010 



and w 



-0.93lo']^5. The la, 2a and 3a contours of Qko and w are shown in Fig. O By 



0.277+°:°i and 



fitting the ACDM model to the combined data, we get x = 333.9, Qmo 
fifco = 0.0002 ± 0.0081. The la, 2a and 3a contours of Qmo and Qko are shown in Fig. [71 
By fitting the DGP model to the combined data, we get = 345.7, Qmo 
Qko = 0.018 ± 0.009. The la, 2a and 3a contours of Qko and w are shown in Fig. [HI 



0.268l|];[J^^ and 



Next, if we use constraints on Q^no and flko from pre vious step and fit + 



to the growth factor data f{z) alone from 
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58 
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61 



01, 



(7-4/7)0^ 



6^ 



64( 1 . we can 



get a constraint on the growth index 7. For the curved ACDM model with the best fit 
value Q^o = 0.277 and Qko = 0.0002, we find that x^ = 4.54 and 7a = O.QbtoH- The 
theoretical value 700 = 6/11 = 0.55 is consistent with the observation at the la level. For 
the DGP model with the best fit value Qmo = 0.268 and Qko = 0.018, we find = 5.79 and 
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iDGP = 0.53lo;i2 which is not consistent with the theoretical value 700 = 11/16 = 0.6875 at 
the la level. 

We also fit Vt1 + {-f- 4/7)^^ to the combined Type la, the BAO, the CMB, the Hubble 
parameter H{z), and the growth factor f{z) data, and allowing for the three parameters 
to vary. We obtain the following constraints on the parameters Qmo, ^ko and 7. For the 
curved ACDM model, we find that = 338.4, n^o = 0.276troliWtomlW^ ^ko = 
-0.0002l°:°[!^^(la)t°:°J?^^2a), and 7 = 0.64l°:^^(la)+°:I^(2a). For the DGP model, we find 
that = 351.47, n^o = 0.2Q8t'of2li^cT)t'o'olli^cx), il^o = 0.017^no9(l^)^o:o}3(2^), and 
7 = 0.52to;27(lc")to;35(2cr). Proceeding in this second way shows that current data that we 
used is not enough in order to put conclusive constraints on 7 and the underlying theory of 
gravity; however, this is likely to change with future data from future high precision missions 
such as PLANCK, JDEM, LSST, and others. 



V. DISCUSSIONS 



We find that the simple analytical expression fa = fij^ + (7 — A/7)Qk can be used to 
approximate very well the growth factor / when the curvature is not zero and this is the 
main result of the paper. The curvature parameter is included, as shown above, into the 
function fa which is found to provide a better fit to the growth factor than Q'^ (with 
curvature included in 7) when compared to the growth factor evaluated numerically from 
the differential equation (this is shown in Figs 3-5). We also find that, with fa as given 
above, the asymptotic value of 7 is the same as that in flat space and it provides distinctive 
information about a dark energy model versus a modification of gravity. For example, 
7 ^ 0.55 for the ACDM model and 7 = 0.6875 for the DGP model. 

We explored then some constraints on 7 using the form of fa above and current data 
from Type la SNe, BAO, WMAP5, H{z) data, and growth factor data. Unlike previous 
analyses which used the shift parameter for the WMAP5 data and the A parameter for the 
BAO data [491], here we used the shift parameter the acoustic scale I a, the redshift 



54| . and the parameter rs{zd)/ Dv{z) at 



and their covariance matrix for the WMAP5 data 
two different redshifts for the BAO data. 

First, we fit the models to the combined Type la SNe, BAO, WMAP5, and H{z) data 

h0.025 Q 
-0.023' '^'fco 



and obtained fi^o = 0.277:^[^:|]23, fi^o = 0.0002±0.0081 and x^ = 333.9 for the curved ACDM 
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model, and n^o = 0.268l[^:S^^, Qko = 0.018 ±0.009 and = 345.7 for the DGP model. The 



DG 



49 



mode^ 
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77 



is found strongly disfavored by the data, a result consistent with other analysis 
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SOj. Comparing with the WMAP5 result |54|, we see that the error bar 



becomes a little bigger with the addition of the H{z) data. The reason is due to the larger 
uncertainties in the H{z) data. Now, using the best fit values for Q^o and Qko from the 
previous step and the growth data of f{z), we obtained constraints on the growth index 
7. For the curved ACDM model, we find that 7a = 0.65lo[i5, which is consistent with the 



38 



40 



491. For the DGP 



theoretical value 0.55. This result is also consistent with those in 
model, we find that ■Jdgp = ^■^'^-oaI- The theoretical value 7 = 0.6875 lies outside the la 
bounds thus again disfavoring the DGP model. 

Next, we fit fa to all the data sets by allowing the three parameters Qmo, and 7 to vary. 
For the curved ACDM model, we find that = 338.4, Q^n = 0.276l^-^^4(la)+°°^^' 



fcO 



-0.0002 



+0.0095/ 
-0.00861 



la)toful{2a), and 7 = OMtoMW^oAtW. For the curved 



'mO 

1+0.48/1 \+0.73/ 



-0.024 v-^^; -0.033 



DGP model, we find that = 351.47, n^o = 0.268I^:^^^(la)I^;^^;^(2cr), fifco = 
OM7tofol{la)1:ofiU'^a), and 7 = 0.52l°:^°(la)t|]:^^(2(T). Thus, the current data is still 
not sufficient for us to obtain conclusive constraints, when we fit all these three parameters 
simultaneously. But this is likely to be improved substantially with future high precision 
missions. 

Finally, the inclusion of the curvature in future analyses of cosmological data in order 
to constrain the origin of cosmic acceleration was recently discussed and stressed in the 
literature, e.g. Q, 81|- This includes the question of distinguishing dark energy from 



modified gravity using the growth data in the presence of spatial curvature. We find here 
that the analytical expression fa = ^l'^ + (7 ~ 4/7)fifc approximates very well the growth 
factor / when curvature is added and should be useful for such future analyses. 
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